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ABSTRACT 

We present Keck/OSIRIS 1.1-1.8 /jm adaptive optics integral field spectroscopy of the planetary- 
mass companion to GSC 06214-00210, a member of the ^5 Myr Upper Scorpius OB association. We 
infer a spectral type of L0±1, and our spectrum exhibits multiple signs of youth. The most notable 
feature is exceptionally strong Pa/3 emission (i?VL=-11.4±0.3 A) which signals the presence of a 
circumplanetary accretion disk. The luminosity of GSC 06214-00210 b combined with its age yields 
a model-dependent mass of 14±2 Mj up , making it the lowest-mass companion to show evidence of a 
disk. With a projected separation of 320 AU, the formation of GSC 06214-00210 b and other very low- 
mass companions on similarly wide orbits is unclear. One proposed mechanism is formation at close 
separations followed by planet-planet scattering to much larger orbits. Since that scenario involves a 
close encounter with another massive body, which is probably destructive to circumplanetary disks, 
it is unlikely that GSC 06214-00210 b underwent a scattering event in the past. This implies that 
planet-planet scattering is not solely responsible for the population of gas giants on wide orbits. More 
generally, the identification of disks around young planetary companions on wide orbits offers a novel 
method to constrain the formation pathway of these objects, which is otherwise notoriously difficult 
to do for individual systems. We also refine the spectral type of the primary from Ml to K7 and 
detect a mild (2-cr) excess at 22 /j,m using WISE photometry. 

Subject headings: stars: pre-main sequence — stars: individual (GSC 06214-00210) - planetary sys- 
tems: formation 



1. INTRODUCTION 

Direct imaging with adaptive optics on large telescopes 
is beginning to reveal the orbital architecture and de- 
mographics of extrasolar planetary systems on wide or- 
bits (>10 AU). Several populations of planetary- mass 
companions are emerging in this nascent field, includ- 
ing giant planets residing in debris disks at moderate 
separ a tions of <100 AU (HR 8799 bede: IMarois et al" 
20081 IMarois et all 120101: Foma lhaut b: IKalas et a! 
20081: B Pic b: lLagrange et all I2009D and planetary 



mass objects on extremely wide orbits of severa l hun- 
dred AU (e.g., C HXR 73 B: ILuhman et all [20061: 
1RXS J1609-2105 b: ILafreniere et al.H 2008t Ross 458 C: 
IGoldman et al.ll20ialsdioldl2010n . 

Despite their growing numbers, the formation mecha- 
nisms of these companions remain obscure. Three plau- 
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sible (and non-mutually exclusive) routes have been pro- 
posed : core accretion plu s gas capture dPollack et al 
1996] lAlibert et al.l 12005ft . disk instability (jCameron 
1978; Boss! I1997D , an d direct co llapse from molecular 
cloud fragmentation (jBatel 12009ft . In situ formation 
through core accretion is unlikely for companions on wide 
orbits because the timescale to grow massive cores at 
these separations is longer than the o bserved lifetimes 
of pr o toplanetary disks (~5 Myr; e.g ., Hernande z et al.l 
[20071: lEvans et al.l 12001 ICurrie et al.l 120091) . On the 
other hand, models of disk instability have succeeded 
in forming gas giants between ~20-100 AU, with spe- 
cific results depending on the initial c onditions and 
physical assumptions of the simulations ( Rafikoyj 120051 : 
IStamatellos k WhitwortIdl2lM lBossll2011ft . 

There is debate about whether disk instability can 
account for planetary-mass co mpanions a t separa- 
tions of several hundred AU (IBossI 120061: iRafikovl 
20071: iNero fc BiorkmaiJl200fl: IStamatellos k Whitworthl 
200a iDodson-Robinson et aLlhoOft" iKratter et al l 120101: 
iRice et alj2010UBolev et al.l2010HBaruteau et al.ll2011ft . 
Some of the conflicting results arise from different ap- 
proaches to modeling disk heating and cooling as well as 
uncertainties in the initial disk masses and surface den- 
sities, which are poorly constrained by observations dur- 
ing the embedded Class and I protostellar phases when 
this mechanism is most likely to occur. Additionally, 
environmental factors such as envelope accretion onto 
the disk (both steady-state and epis odic) are only be- 
ginning to be included in sim ulations (jVorobvov k Basul 
l2010t IStamate llos et al.|[2011l ). Despite these difficulties, 
there is some evidence that disks around Class objects 
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can be both massi ve (>0.2 M w ) and large (>200 AU; 
Eisner et al.l 120051: Hirgensen et all 120051; lEnoch et all 
2OO9HJ0rgensen et alJ l2009) , leaving open the possibility 
that disk instability can form planets at wide separations. 

There have been fewer theoretical studies focusing on 
the formation of gas giants by direct collapse from molec- 
ular clouds. In this scenario, planetary-mass objects 
form as ejec ted by-prod ucts of a fragmenting pre-stellar 
cloud core. iBatd ()2009f ) used numerical simulations to 
follo w the gravitationa l coll apse of a molecul ar cloud (see 
also iBate et al.l 120021 and IBate et al.l I2003Q and found 
that all stellar and substellar objects begin as opacity- 
limited fragments with masses of a few times that of 
Jupiter and subsequently accrete gas, increasing their 
mass over time. In these simulations low-mass brown 
dwarfs and planetary-mass objects form directly from 
dense filamentary cloud structures and from instabilities 
in disks; in both cases they begin as dynamically unsta- 
ble multiple systems and are ejected from dense regions 
of gas, halting strong accretion and limiting their masses 
below the hydrogen burning minimum mass. The forma- 
tion of gas giants at extreme separations may have also 
proceeded in this fashion, perhaps by being ejected to 
wider orbits on faster timescales than brown dwarfs be- 
fore appr eciable accretion h as occurred. Alternatively, as 
noted by IBate et al.l (|2003). preferential accretion from 
the primary could also explain the population of both 
brown dwarfs and planetary-mass companions on orbits 
of several hundred AU. In this framework of fragmenta- 
tion plus ejection, planetary-mass objects on wide orbits 
represent the low-mass tail of brown dwarf companion 
formation. 

In addition to these in situ formation models, other 
explanations for giant planets on wide orbits involve 
formation at close separations (perhaps through core 
accretion or disk instability) and subsequent orbital 
evolution to large separations. Several possibilities 
have been proposed, including outward scattering from 
dynamical interactions with another massive planet 
(IBqssI 120061: iDebes fc Sigurdssonl l200l iScharf fc Menoul 
120091 : iVeras et al.l 120091 ) and outward resonant migra- 
tion with another planet while still embedded in a disk 
(|Crida et al J 1200 9). These myriad possibilities make it 
difficult to identify the formation mechanisms of individ- 
ual systems discovered by direct imaging. 

The value of this population extends beyond inform- 
ing formation scenarios; the direct detection of thermal 
photons enables detailed studies of their atmospheres. In 
particular, spectroscopy of young low-mass companions 
is providing insight into the influence of surface grav- 
ity on th e atmospheric propertie s of lo w-temperature 
objects. iMetchev fc Hillenbrand! (j2006l ) find that the 
young brown dwarf companion HD 203030 B has an 
earlier spectral type than expected from its evolution- 
ary model-derived temperature, suggesting that grav- 
ity may impact the transition from L-type dwarfs to 
T-type dwarfs. More recent photo metry and spec- 
trosc opy of the young (^30 Myr; iZuckerman et al.l 
120 lit ) HR 8799 planets indicate that they have unusu- 
ally thick photospheric clouds, possibly signaling that 
thick clouds is a general phenomenon of young gas gi- 
ants dMarois et alj200l iBowler et al.ll20iq iCurrie et al.l 
120111: iMadhusudhan et all 120111: IBarman et al.l l2011al ). 
The emergent spectrum of the ~5 Myr planetary-mass 



compa nion 2M1207b shows similar signs of a dust y atmo- 
sphere (jSkemer et al.ll2011h IBarman et al.ll2011bO . There 
is a a growing need to explore this unexpected correla- 
tion with more objects spanning a range of gravities and 
temperatures (or equivalently, masses and ages). 

Here we present near-infrared spectroscopy of the re- 
cently dis covered planetary-mass companion G SC 06214- 
00210 b (jKraus et al.l [20081 ; llreland et afll20lTI ). which 
orbits a member of the Upper Scorpius OB association 
(^5 Myr) at a projected separation of 320 AU. In §2 
we present J- and iJ-band spectroscopy of the compan- 
ion as well as optical and near-infrared spectroscopy of 
the primary. We describe the spectroscopic properties 
of GSC 06214-00210 b, atmospheric model fits, predic- 
tions from evolutionary models, and analysis of the pri- 
mary star in §3. In §4 we examine the validity of planet- 
planet scattering as a migration scenario for GSC 06214- 
00210 b. Finally, §5 provides the conclusions of our work. 

2. OBSERVATIONS 

2.1. Keck/OSIRIS J- and H-Band Spectroscopy of 
GSC 06214-00210 b 

We observed GSC 06214-00210 b on 9 July 2010 
with the OH-Suppressing Inf rared Imaging Spectrograph 
(OSIRIS: Earkin et al.l l2006) integral field unit with nat- 
ural guide star adaptive optics on Keck II. The weather 
was photometric with good seeing (0'.'4-0'.'6 according to 
the DIMM on CFHT). We observed GSC 06214-00210 b 
in the Jbb (1.180-1.416 ,um) and Hbb (1.473-1.803 /im) 
bandpasses with the 0'/05 pix" 1 plate scale, resulting in 
a lenslet geometry of 16x64, a field of view of 0'/8x3'.'2, 
and a resolving power (R = X/SX) of ^3800. The sep- 
aration of GSC 06214-00210 b from the primary (2'.'2) 
is large compared to the measured FWHM of the com- 
panion (O'/ll in Jbb; 0'.'13 in Hbb) and the seeing disk 
(~0'/5), so contamination from the primary is negligible. 
We obtained a total of 40 min in Jbb (5 min/exposure 
x 4 nodded pairs) and 30 min in Hbb (5 min/exposure 
x 3 nodded pairs) with 1" offsets along the long axis 
of the field of view. The airmass ranged from 1.35- 
1.64 during the observations. Immediately following our 
science observations we targeted the nearby A0V star 
HD 148968 at an airmass of 1.50-1.67. We also acquired 
sky frames in both filters before GSC 06214-00210 and 
after HD 148968. 

The raw data were flat fielded, sky subtracted, cleaned 
for cosmic rays and bad pixels, assembled into 3D data 
cubes using the appropriate rectification matrices, and 
wavelength calibrated using version 2.3 of the OSIRIS 
Data Reduction Pipeline. Figure Q] shows an example 
of a collapsed J-band data cube. The science target 
and standard star spectra were extracted with aperture 
radii of 3 and 4 spaxels, respectively. The individual 
spectra from each band were then scaled to a common 
level and median combined. Spectral measurement un- 
certainties were determined by computing the standard 
errors about the median. The spectra were telluric cor- 
rected with the xtellcor_general ro utine in the Spex- 
tool spectroscopic r eduction package (jVacca et al.l 120031 : 
iCushing et al .1120041 ) . We tested sky subtraction with sky 
frames (A-sky) and with the nodded science data (A-B) 
by reducing the data both ways; the influence on the fi- 
nal spectrum was minor. We chose the former (A-sky) 
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for our final spectrum. 

The J- and -ff-band spectra were independently flux 
calibrated us i ng flu x ratios for GSC 06214-00210 from 
llreland et "all (120111) together with 2 MASS photometry o f 
the primary (jSkrutskie et al.ll2006D . llreland et all ([20111 ) 
report flux ratios in the MKO system, so we use the 
2MASS-MKO conversions from iLeggett et all (|2006l ) to 
arrive at MKO apparent magnitudes for GSC 06214- 
00210 b of J=16.25±0.04 mag and 7i=15.55±0.04 mag. 

2.2. IRTF/SpeX Near-Infrared Spectroscopy of the 
Primary GSC 06214-00210 

We observed the prima ry star GSC 06214-0 0210 with 
the SpeX spectrograph (|Ravner et al.l [2003) in short 
wavelength cross dispersed mode (SXD) at the NASA 
Infrared Telescope Facility (IRTF) on 29 April 2011 UT. 
The seeing reported by the CFHT DIMM was between 
0'.'4-0'. / 6, and there were light cirrus clouds during the 
observations. We used the C/3 slit aligned to the par- 
allactic angle, yielding a resolving power of —2000, and 
obtained a total of 4 min of data (30 s/exposure x 4 
nodded pairs) by nodding along the slit in an ABBA 
pattern. After our science observations we targeted the 
nearby A0V star HD 144925 and obtained calibration 
frames at a similar air mass. The data were reduced with 
versio n 3.4 of Spextool ([Vacca et al.ll2005 iCushing et al.l 
l200l . 

We also acquired low-resolution (i?— 100) spectroscopy 
of the primary with IRTF/SpeX in prism mode on 12 
May 2011 UT. The weather was poor with complete cir- 
rus coverage so we used the opportunity to verify the 
spectral slope of our SXD spectrum. We used the 0'.'8 
slit and obtained 160 s of data (20 s/exposure x 4 nod- 
ded pairs) at an airmass of 1.38. Immediately after we 
observed HD 144925 for telluric calibration. 

2.3. UH 2.2m/SNIFS Optical Spectroscopy of 
GSC 06214-00210 

We obtained an optical spectrum of GSC 06214-00210 
on 16 May 2011 U T under clear con ditions using the 
SNIFS instrument ([Lantz et al.l 12004ft on the UH 2.2m 
telescope. SNIFS is an optical integral field spectro- 
graph with i?~ 1000-1300 that splits the signal with 
a dichroic mirror into blue (-3000-5200 A) and red 
(—5200-9500 A) channels. The images are resamplcd 
with microlens arrays, dispersed with grisms, and focused 
onto blue- and red-sensitive CCDs. 

A single 100 s exposure of the science target was suf- 
ficient to achieve high S/N (—200). The processing 
was performed with the SNI FS data reduction pipe line, 
which is described in detail in lAldering et al.l ([20061 ) and 
IScalzo et al.l §010). The processing includes dark, bias, 
and fiat-field corrections; assembling the data into red 
and blue 3D data cubes; cleaning them for cosmic rays 
and bad pixels; sky subtraction; extracting the spec- 
tra using a semi-analytic PSF model; and wavelength- 
calibrating the spectra with arc lamp exposures taken at 
the same telescope pointing as the science data. Cor- 
rections for instrument response, airmass, and telluric 
lines are based on observations of the Feig e 66 st andard 
star using calibrated observations in lOkd ([1990). The 
far ends of the blue and red channels have low QE so to 
avoid these regions we trim our final spectrum to 3300- 



4900 A and 5170-8700 A. 

3. RESULTS 

3.1. Properties of GSC 06214-00210 b 

3.1.1. Spectral Properties and Classification of 
GSC 06214-00210 b 

Our OSIRIS spectrum of GSC 06214-00210 b is pre- 
sented in Figure The most striking feature in the 
data is the strong Pa/3 emission (EW "=-11 .4±0 .3 A) at 
1.282 /im. The feature is present is each pair of dithered 
J-band observations of GSC 06214-00210 b and spans -9 
spectral elements (Figure [21 inset), so it is not a result 
of a cosmic ray or bad pixel. One of our J-band OSIRIS 
data sets include the primary (seen in Figure[TJ the rest of 
data are dithered to avoid the star) and no Pa/3 emission 
is observed in that spectrum of the star. Pa/3 emission is 
also absent in our medium resolution SXD spectrum of 
the primary obtained —8.5 months later. We therefore 
rule out the possibility that the emission observed in the 
companion is a result of contamination from the primary. 

Pa/3 emis sion in young stars can arise from accretio n 
or outflows (|Folha fc Emersonl200lllWhelan et al.l2004ft . 
both of which imply the presence of a circumplanetary 
disk around GSC 06214-00210 b. The presence o f a dis k 
is bolstered by photometry from llreland et al.l (|2011ft . 
who found a red K-L' color compared to field M and L 
dwarfs and suggest the excess may originate from ther- 
mal disk emission. Although other parameters such as 
metallicity, surface gravity, and dust can also affect the 
K-L' color of late-type objects, these are unlikely to be 
the origin of the red color for GSC 06214-00210 b based 
on the expected near solar-metallicity of USco mem- 
bers, a comparison to colors of field objects and giants 
([Ireland et al.ll20TlT l . and the minor influence of dust on 
K — L' colors at the M/L transition predi cted by models 
(|Baraffe et al.ll2003l: IChabrier et al.ll2000ft . We note that 
the K-L' color of 1.2 mag is redder than many other 
youn g brown dwarfs with excesses attr ibuted to disks 
(e.g.. lLiu et"all[200l lAllers et al.l [2001 ). W e desc ribe 
the inferred accretion rate and luminosity in £13.1.41 and 
we discuss the implications of a disk in the context of the 
formation of GSC 06214-00210 b in gU 

Other notable features in our spectrum include absorp- 
tion at -1.20 ^m from FeH, 1.244 and 1.253 /im K I 
lines, deep H2O steam bands at —1.4 and —1.8 /im, and 
FeH bandheads in the H band near 1.6 /im (Figure [5]). 
All of these point to a late-M or early-L spectral type 
(|Cushing et al.ll2005h . 

Further classification of GSC 06214-00210 b warrants 
some care since there is no universally adopted near- 
infrared classification scheme for young late-type stars 
and brown dwarfs. The classification of young brown 
dwarfs in the literature is inconsistent and often a mix of 
qualitative (co mparative) and quantitative (index-based) 
systems (e.g.. [Reid et al.l Il9 95; Kir kpatrick et all 11999; 



Martin et all 1999t iReid et al.ll200H iGeballe et all 
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iRiddick et alJl2007ft ; when multiple systems are used the 
adopted spectral type is usually a (subjectively weighted) 
average from several schemes. 

We classify GSC 06214-00210 b primarily based on 
comparisons to near-infrared spectra of field objects and 
late-type free-floating members of Upper Scorpius (Fig- 
ure [3]). The advantage of using field objects is that they 
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are optically typed and the advantage of using USco 
members is that they have the same age as GSC 06214- 
00210 b. Spectra of field objects origi nate from the IRTF 
Spectral Libr ary dRavner et al.ll2009( ) and those of USco 
objects from lLodieu et al l (|2008t ). Many of the USco 
spectra have modest S/N so we first Gaussian smoothed 
them from their native resolution of i?^1700 to i?~1000 
and cleaned them for strong single pixel outliers proba- 
bly originating form imperfect bad pixel or cosmic ray 
removal in the original reduction. For the comparative 
analysis in FigureGjjthe OSIRIS spectrum of GSC 06214- 
00210 b was smoothed to the appropriate resolving power 
for each sample. 

The K I lines in GSC 06214-00210 b are weaker than 
those in M8-L2 field objects, which i s a common sig- 
natur e of youth (IGorloya et all 120031 : IMcGovern et al.l 
12004 IKirkpatrick et al.l I2006TT When normalized to 
1.29-1.32 /im, as in Figure G2 the best matches to the 
J band are M8-L0.5 dwarfs. The J-band in GSC 06214- 
00210 b is significantly bluer compared to the field LI 
and L2 dwarfs. Together the joint J+H spectra appear 
to best fit the L0.5 dwarf, although the 1.6-1.65 /im FeH 
features are somewhat stronger in the field object. 

The best USco matches are M9-L2 objects in the J 
band and M9-L1 objects in the H band. Interestingly, 
the 1.6 /im FeH features appear to be slightly stronger 
and the overall H band less triangular in GSC 06214- 
00210 b compared to the USco members. The depth 
of the K I lines are similar for GSC 06214-00210 b 
and the M8-L2 sequen ce. Note that the sp ectral types 
of some objects from lLodieu et al.l (|2008|) appear to 
be 1-3 subtypes later than their opt ical classifications 
dHerczeg et al.ll2009t iBiller et al.ll201lD . 

lAllers et al.l (|2007[ ) define a spectral index from 1.50- 
1.57 /im based on the depth of the H 2 absorption that 
correlates well with optical spectral types from M5-L0 
and is independent of gravity. This index yields a spec- 
tral type of M9.5t^ for GSC 06214-00210 b, which is 
consistent with the comparison to field objects and USco 
members. The errors incorporate spectral measurement 
uncertainties and intrinsic scatter in the relation. Alto- 
gether we assign GSC 06214-00210 b a spectral type of 
L0±1. No te that this agrees well with the estimate of 
M8-L4 bv llreland et al.l (|201ll ) based on colors. 

3.1.2. Effective Temperature and Bolometric Luminosity 

We fit the solar metall icity BT-Settl grid of model at- 
mospheres (2010 version; lAllard et"aLll2010D to our spec- 
trum of GSC 06214-00210 b. The BT-Settl models in- 
corporate updated water opacity sources, revised solar 
abundances, and a new advanced treatment of dust for- 
mation using radiation hydrodynamic simulations. This 
new grid reproduces the observed colors and SEDs of 
late-M and L dwarfs b etter than the Ame s-Cond and 
Ames-Dusty models of lAllard et al.l ([2001D . which ex- 
amined the limiting cases of atmospheric dust forma- 
tion. The grid spans effective temperatures between 
2000-3000 K (AT cff =100 K) and gravities between 2.5- 
5.5 dex (cgs; Alog g=0.5). The fitting procedure reli es on 
Y 2 minimizat i on as detailed in lCushing et al.l ()2008l ) and 
IBowler et all (|2009). To improve the S/N of the data we 
Gaussian smoothed our OSIRIS spectrum from its native 
resolving power of i?~3800 to i?~2000 and propagated 



the measurement errors through the convolution. We 
smoothed the models to i?~2000 and fit them to the data 
in a Monte Carlo fashion by randomly generating syn- 
thetic observations based both on the spectral measure- 
ment uncertainties and the J- and -ff-band photometric 
uncertainties ( EI2.1j) used for flux calibration. For each 
trial we determine the best-fitting spectrum and save the 
X 2 values from fitting every model in the grid. Figure [4] 
displays \ 2 contour maps averaged over 100 Monte Carlo 
trials. Since the spectra are flux calibra ted and the dis- 
tance (d) is constrained t o 145±14 pc (|de Zeeuw et al.l 
[19991: llreland et afll20Tl . the radius (R) can be cal- 
culated from the multiplicative fac tor (R 2 /d 2 ) used t o 
scale the models to the data (see IBowler et al.l [2009) . 
Figure 2] shows the contour plot of the inferred radius 
for GSC 06214-00210 b, representing averages from the 
Monte Carlo trials and incorporating the uncertainty in 
the distance estimate. 

The J (1.17-1.36 //m). H (1-45-1.83 /im), and J+H 
(1.17-1.83 /mi) spectral segments were fit separately. To 
avoid the Pa/3 emission line the 1.280-1.285 /im region 
was excluded from the fits. The best-fitting model to the 
J band is T off =2700 K/log g=4.5 with #=1.22 i?j up 0] 
The models match the depth of the K I lines and 
1.35 /im H2O depth reasonably well, but they fail to 
reproduce the 1.2 /im FeH feature and the continuum 
level from 1.27-1.31 /im. In the H band the models 
provide a better match the data. The best-fit model 
is T Gff =2700 K/log g=4.0 with #=1.36 # Jup , which is 
similar to the result from the J band. 

A fit to the combined J+H region produces a different 
result: T eff =3000 K/log g=2.5 and #=1.02 i? Jup . The 
quality of the fit is rather poor as the J band continuum 
slope deviates dramatically from the data and the depth 
of the K I lines are poorly reproduced. The H band, 
however, is accurately reproduced by the model. In this 
case the best-fitting model is at the edge of the model 
grid so it may not represent the global minimum. The 
2700 K/4.0 dex model that matches the individual bands 
overestimates the J-band flux and underestimates the H- 
band flux when fit to the entire spectrum. 

Another approach to infer effective temperature is 
to use predictions from evolutionary models based 
on an object's age and bolometric luminosity. The 
age of the Upper Scorpius association has long been 
known to be ^5 Myr, primarily constrained from fitting 
isochrones to known members in the HR diagr am (see 
iPreibisch fc Mamaiekll200l iSlesnick et al.ll2008D . Infer- 
ring the age of a cluster from the HR diagram can be 
problematic, however, because pre-main sequence evo- 
lutionary models may have significant systematic errors 
(|Hillenbrand fc White! [20041: iKraus fe Hillenbrand! 1200911 
perhaps caused in part by the influence of episodic ac - 
cretion ([Baraffe et all 120091: iBaraffe fc Chabrierl 1201(1 . 
With these uncertainties in mind, we adopt an age of 
5±2 Myr for GSC 06214-00210 b. 

We calculated the bolometric luminosity by integrat- 
ing an artificial spectrum constructed from the flux cali- 
brated OSIRIS J and H band spectral segments com- 
bined with a model spectrum for wavelengths shorter 
than the J band segment (0.001-1.15 /im), between the 

7 The standard value for Jupiter's radius is 71,492 km 
ILindal ct al. 1981), which is the equatorial radius at 1 bar. 
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J and H segments (1.36-1.45 (an), and longer than the 
H band segment (1.83-1000 //m). At each model-data 
interface the model was scaled to the data to create a 
continuous spectrum. Uncertainties in the spectra and 
flux calibration were accounted for in a Monte Carlo fash- 
ion. To test the sensitivity of the resulting luminosity on 
the input atmospheric model temperature we calculated 
the luminosity using three BT-Settl models with temper- 
atures of 2000 K, 2500 K, and 3000 K, and a gravity of 
4.0 dex. The influence on the resulting luminosity was 
negligible (~0.01 dex). Based on 1000 Monte Carlo trials 
we derived a bolometric luminosity of -3.1±0.1 dex for 
CSC 06214-00210 b. 

To further verify our luminosity calculation we use 
the fsT-band bolometric correcti o n cali brated with field 
objects from iGolimowski et ahl (|2004| ). We arrive at 
a luminosity of -3.04±0.10 dex, where the uncertainty 
takes into account errors in the distance, spectral type, 
photometry, and intrinsic scatter in the relation. In 
addition we compute the luminosity using parame- 
ters predicted from evolutionary models (see below; 
T c fr=2200 K/log g=4.0). We flux calibrated the BT- 
Settl, Ames-Dusty, and Ames-Cond (|Allard et al.ll2001f ) 
models to the mean of the J-, H-, and if-band flux cal- 
ibration scaling factors, accounting for photometric un- 
certainties in a Monte Carlo fashion. The resulting lumi- 
nosities are nearly identical at -3.11 ±0.09 dex, only dif- 
fering by —0.01 dex among the models. The radii inferred 
from the scaling factor and the distance (see above) are 
also virtually the same at 1.89±0.18 i?j up . 

Figure [5] (top panels) displays interpolated temper- 
ature tracks as a function of luminosity and age for 
various evolutionary models. The left panel shows the 
Lyo n Cond and Dusty m odels of iBaraffe et al.l (|2003f ) 
and iChabrier et ahl (|2000l ). which demonstrate the lim- 
iting effects of photospheric dust forma tion (forma- 
tion and settling vs. complete retention; lAllard et al.l 
l2001[h The right panel show s the evolutionary models 
of ISaumon fc Marlevl (|2008f) for photospheres with no 
clouds ( "nc" ) and those with significant amounts of dust 

("/sod = 2"). 

We determined effective temperatures in a Monte 
Carlo fashion by interpolating the evolutionary models. 
The Lyon models begin at 1 Myr and the Saumon & Mar- 
ley models begin at 3 Myr, so we excluded ages younger 
than those in the analysis. (The resulting PDFs 
for the age are truncated Gaussians in linear space). 
The Cond and Dusty models yielded nearly identi- 
cal results of 2160±85 K for 10 4 Monte Carlo draws. 
The Saumon & Marley models were similar, yielding 
2200±100 K for the case with no clouds and 2185±120 K 
for the cloudy version. These values are substantially 
lower than the those of the best-fitting model atmo- 
spheres. 

The temperatures inferred from the atmospheric mod- 
els (—2700-3000 K) are warmer than previous determi- 
nations of young M/L transition objects in the literature 
by -300-500 K (e.g ., IBeiar et al.|[200l IKuzuhara et~aTl 
20111: iLuhman I2004T) . A similar result was obtained by 



els, which are thought to be more reliable because they 
are less sensitive to m issing or incomplete opacity sources 
(|Chabrier et al.l 120001 ) . Likewise, the radii inferred from 
the atmospheric model fitting arc systematically smaller 
than the value of — 1.95 i?j up predicted by evolutionary 
models f i)3.1.4p . which is probably a result of overesti- 
mated effec t ive te mperatures in the fits. 

ILuhman! (|1999[ ) developed an effective temperature- 
spectral type scale for young M dwarfs intermediate 
between those of dwarfs and giants. The calibration 
was defined to ensure that components of the quadru- 
ple system GG Tau and members of the young cluster 
IC 348 were aligne d on the same isochrone from the 
IBaraffe et al.l (119981) evoluti onary models in the HR dia- 
gram. ILuhman* et al.l (120031 ) revised the scale for M8 and 
M9 spectral types based on the latest-type members of 
IC 348 and Taurus. According to Luhman et al., the ef- 
fective temperature of a young M9 object is —2400 K. To 
estimate temperat ures beyond the Luhman et al. scale, 
lAllers et al.l (|2007f ) subtract relative offsets taken from 
the temperature scale of field objects. This amounts to 
-92± 175 K for field M9-L0 ob jects based on the rela- 
tion of IGolimowski et al.l (|2004T > (the uncertainty is from 
the rms scatter of the SpT-T e g relation). Assuming ap- 
proximate uncertainties of 100 K for the Luhman et al. 
scale, this yields a temperature of — 2310±200 K, which 
is consistent with the evolutionary model predictions but 
disagree with those from the atmospheric model fitting. 

3.1.3. Mass 

The mass of GSC 06214-00210 b wa s first estimated 
to be -12 M Jup bv IKraus et~aTI (|2008D from M K mea- 
surements. It was r ecently updated to — 12-15 Mj up by 
llreland et al.l (| 2 1 If) based on JHK colors and evolu- 
tionary models. We refine the mass determination us- 
ing our luminosity measurement and a variety of evolu- 
tionary models (Figure [5l bottom panels). We use the 
same procedure to determine the mass from evolutionary 
models as we use in i)3.1.2l to determine temperature. 
The Cond and Dusty models yield 13.6±2.4 Mj up and 
14.2±2.4 Mj up , and the Saumon & Marley cloudless and 
cloudy models yield 14.1±1.9 Af Jup and 14.4±1.8 M Jup . 
These values are in close agreement with previous esti- 
mates and hug the border of the brown dwarf/planetary- 
mass limit as defined by the deuterium-burning limit 



ass limit as denne d by tne c 
T3 M T „ r : ISpiegel et al.ll20Tl . 



Dupuv et aiT (j2010) 

in their analysis of resolved late- 
M binaries in the field with dynamical masses. They 
found that the most widely used atmospheric models sys- 
tematically yield effective temperatures —250 K higher 
than the temperatures predicted by evolutionary mod- 



3.1.4. Accretion 

In a study of accretion diagnostics in young brown 
dwarfs. iNatta et all (|2004f ) found that Pa/3 luminosity is 
well correlated with accretion luminosity at low masses. 
We use their empirical relation to derive the accretion lu- 
minosity and mass accretion rate for GSC 06214-00210 b. 

The equivalent width of the Pa/3 emission line is - 
11.4±0.3 A, and the flux from that line is 1.12±0.03 x 
10~ 18 W mT 2 . Assuming a distance of 145±14 pc we 
find log(L Pa/3 /L )= -6.14 ± 0.08 for the Pa/3 line lumi- 
nosity. Applying the empirical relationship between Pa/3 
lumin osity and accretion luminosity from INatta et al.l 
(|2004 Equation 2) yields an accretion luminosity of 
log(L acc / Lq)= —4.4 ± 1.3. The mass accretion rate 

(M) and the accretion luminosity are related through 
M=L acc R/GM, where R is object's radius and M is 
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its mass. We calculate a radius from evolutionary mod- 
els in the same fashion as in £13. 1 .21 and £|3. 1 .31 and find 
i?=0.20±0.01 Rq. This yields a mass accretion rate of 
log(M)=-10.7±1.3, where M is in M yr _1 . The un- 
certainty is dominated by scatter in the fitted relation 
from Natta et al. and reduces to 0.14 dex when these 
are ignored. 

Mass accretion rates are observed to depend strongly 
on stellar mass and roughly follow an Moi.M'j empiri- 
cal relationship (e. g.. IMuzerolle et al.|l2005b iNatta et al.1 
120061 although see Idarke fc PringldboM . This trend 
spans two orders of magnitude in mass and six orders of 
magnitude in mass accretion rate, although there are over 
two orders of magnitude of intrinsic scatter in the rela- 
tion. Figure E] shows the position of GSC 06214-00210 b 
in the logM-logM* diagram relative to members of var- 
ious star forming regions. The dependence of accretion 
rate on mass is clear, and objects with masses <30 Mj up 
(log(M*/M©)<-1.5) appear to have accretion rates be- 
low ~10~ n Mq yr" 1 . GSC 06214-00210 b has a some- 
what higher accretion rate than the lowe st mass brown 
dwarf s from IMuzerolle et al.l (|2005l ) and IHerczeg et al.l 
( 2009) , although it is consistent with the scatter seen at 
higher masses. 

While we have assumed the emission in GSC 06214- 
00210 b originates from energy released during accretion, 
we briefly examin e whether it could res ult from chromo- 
spheric activity. IShort fc Dovld ([1998D modeled active 
M dwarf chromospheres and found that weak Pa/? emis- 
sion is possible at some chromospheric pressures. In Fig- 
ure [7] we show the relative strength of the Pa/3 EW and 
line luminosity to objects from the INatta et al.l (2006) 
sample of young stars. The Pa/3 emission strength for 
GSC 06214-00210 b is comparable to the strongest ac- 
cretors from that sample and is larger than the vast 
majority of comparison stars. We can also examine 
the correlation between Pa/3 emission and the Ha 10% 
line width, which is a widely-adopted accretion indica- 
tor. Amon g the sample of young brown dwarfs from 
INatta et"a l. (2004), four out of the five that show ev- 
idence of accretion based on their Ha 10% widths also 
exhibit Pa/3 emission. No objects with Pa/3 emission were 
identified as non-accretors, which together with the un- 
usually strong Pa/3 EW supports the notion that the 
observed emission in GSC 06214-00210 b is a result of 
accretion. 

We note that GSC 06214-00210 b is also likely to have 
other emission lines. Ha is a good candidate since it 
correlates well with Pa/3 emission in brown dwarfs. Br7 
is observ ed in some, but no t all, accreting objects. For 
example, INatta et al.l ((2004) found that only two out of 
eight brown dwarfs with Pa/3 emission showed Br7 emis- 
sion. However, Br7 emission does seem to correlate with 
the strength of the Pa/3 line, which increases the likeli- 
hood that GSC 06214-00210 b is emitting at that line. 

3.2. Properties of the Primary GSC 06214-00210 

The primary star GSC 06214-002 10 is a weak-lined 
T Tauri member of Upper Sco (see llreland et al.l 120111 
for a review of the literature). Medium r esolution opti- 
cal spectroscopy bv lPreibisch et al] ([1998) revealed weak 
Ha emission (EW=-\.hl A) which can be attributed to 
chromospheric activity. Our 0.8-2.5 fun SXD spectrum 



is presented in Figure [5J no emission lines are present. 
Compared to lat e-K and early-M dw arfs from the IRTF 
spectral library (|Ravner et al.l [2009) , as shown in Fig- 
ure [9] (right panel), the near- infrared spectrum resembles 
late-K objects more closely than early-M objects. We 
verified that our SXD spectrum was properly reduced by 
obtaining a low-resolution spectrum with SpeX in prism 
mode during a night of poor weather conditions. The 
shapes of the spectra are virtually identical (lower left in- 
set of Figure [8]) , verifying our SXD spectrum and calling 
into questio n the nominal s p ectra l type of Ml originally 
assigned bv lPreibisch et al] (j!998l ). 

Overall, the optical sp ectrum i s best matched by that 
of a K7 dwarf from the iPickled (|1998f) spectral library 
from ~3300-8700 A (Figure U left panel) @ The shape 
of the SED and the depth of the TiO absorption bands 
differ substantially from those of the Ml spectrum. We 
also infe r a spectra l type using the spectral indices de- 
fined by iReid et al.l (|1995f ). which measure the depth of 
various molecular absorption bands in the optical. Our 
value for the Ti05 index (0.916) yields a spectral type 
of K5.5. Since many of the optical absorption bands 
are gravity-sensitive, a slightly earlier spectral type is 
not unexpected using this index-based scheme. Simi- 
larly, the CaH2 index (0.908) indicates an effective tem- 
perature of ^4200 K usin g the SpT-T c s relation from 
IWoolf fc Walierste in (2006). All of these diagnostics sug- 
gest a spectral type of late-K rather than early-M, so we 
revise the spectral type of GSC 06214-00210 from Ml to 
K7±0.5. 

We determine the effective temperature of GSC 06214- 
00210 by constructing its SED from 0.6-22 /im (Fig- 
ure I10p and comparing it to the gri d of solar metallicity 
Phoenix-Gaia model atmospheres (Brott & Hauschildt 
2005). The photometry originate from t he Carlsberg 
Merid ian Catalog 14 (CMC 14) for r' band ([Evans et all 
120021 ); the Deep Near-Infrared So uthern Sky Survey (DE - 
NIS) for i, J, and K s bands ([Epchtein et al] 119971) : 
the Two Micro n All Sky Survey (2M ASS) for J, H, 
and K s bands ([Skrutskie et al.l I2006D : and the Wide- 
Field Infrared Sur vey Explorer (WISE ) for 3.4, 4.6, 12, 
and 22 /urn bands ([Wright et alj|2010f). Z ero point flux 
densities are f rom [stoughton et all (12002ft for CM C14, 
iFouaue et all (120001) for DENIS, IRieke et all ((2008ft for 
2MASS. and lWright et all (pOlOt ) for WISE. The models 
arc flux calibrated using the mean scaling factor from the 
2MASS J, H, and Ks bands. We limit the comparison 
to models with log <?=4.0, which is the approximate sur- 
face gravity from evolutionary models for a K7 star at 
5 Myr. The r' , i, and J bands carry the most weight in 
constraining the effective temperature because the SED 
turns over to the Rayleigh- Jeans tail at A>1.6 /.tm (H 
band). The best match is the 4200 K model, although the 
4100 K and 4300 K models provide decent fits. Warmer 
and cooler model temperatures begin to diverge from the 
photometry, so we assign an uncertainty of 150 K for the 
temperature. 

The WISE 22 /im photometry disagrees with the 
model in Figure [IjJ suggesting a slight excess for the 
primary. In Figure [Tl] we compare the position of 
GSC 06214-00210 to field K7 dwarfs and other USco 

8 We use the modification to the MK classification scheme that 
assigns K7 as an intermediate type between K5 and M0. 
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members in the Wl-WA vs J-Kg diagram. To isolate a 
sample of field K7 dwarfs we used a compilatio n of the re- 
vised Hipparcos catalog ([van Leeuwenl 120071 ) generated 
by E. Mamajek (private communication) with spectral 
types and y-band mag nitudes from the original catalog 
(jPerryman et al.l 119971 ). An absolute magnitude cut of 
My > 5 was imposed to exclude giants. This yielded 272 
K7 dwarfs, which we then fed into the WISE Prelimi- 
nary Release Source Catalog query using a search radius 
of 10". (The PSF FWHM for WISE bandpasses range 
between 6-12".) Among the resulting detections, we 
kept only those with the best photometric quality flags 
in all bands (ph_qual="A" or "B"), the best contami- 
nation and confusion flags (cc-flags= "0000"), extended 
source flags consistent with a point source (ext_flg= u 0"), 
and stable variable flags (var_flg<5). This produced 
41 objects, which arc plotted in Figure [TTJ We per- 
formed a similar search of the WI SE Preliminary Cata- 
log for known USco members from Preibi sch k Mamaiekl 
(2008). We queried around 251 USco members with 
GKM spectral types applying the same flags as above; 
67 objects yielded reliable WISE detections. Field 
stars and some USco members have Wl-WA colors of 
— 0.0±0.2 mag, although most USco objects show mod- 
erate to large excess in that color. The W1-W4 color 
of 0.78±0.27 for GSC 06214-00210 indicates a mild (2-cr) 
excess compared to field K7 stars. 

We calculate the luminosity of the primary by integrat- 
ing the flux calibrated 4200 K/log g=4.0 synthetic spec- 
trum, yielding L*=0.38±0.07 L@. The error budget is 
dominated by the uncertainty in the distance (the uncer- 
tainty in effective temperature contributes —0.01 dex). 
In Figure [12] wc show the position of GSC 06214-00210 
in the HR di agram. The pre-ma in sequence evolution- 
ary models of iBaraffe et al.l (|1998| ) for L m i X —H p are dis- 
played from 0.02-1.4 Mq with isochrones from 1 Myr 
to 1 Gyr (black). We also include the L m i x =1.9H p 
1 Myr, 5 Myr, 10 Myr, and 1 Gyr isochrones for higher 
masses (>0.6 Mq). Gray circles show the positions of 
Upper Scorpius low- and int ermediate-mass members in 
iPreibisch k Mamaiekl ([20081 Tables 1 and 2). The po- 
sition of GSC 06214-00210 is consistent with the scatter 
exhibited by other Upper Scorpius members, although 
it sits at a somewhat older isochrone than the nominal 
5 Myr age of the complex. For the prescription in Fig- 
ure[12]with L m i X —H p , the isochronal age is —16 Myr with 
a l-cr range of 13-25 Myr. For L m ; x =1.9.ffp, which is re- 
quired to fit the Sun, the inferred age is — 10 Myr with a 
range of 6-16 Myr. Note that the isochronal age of low- 
mass USco members (—5 Myr) diverges from the cluster 
age based on the L m i X —H p prescription (~10 Myr), and 
instead is more consistent with the age inferred from the 
L m i K =1.9H p grid. 

We estimate the mass of GSC 06214-00210 from 
its position on the HR diagram together with pre- 
dictions from pre-main sequence evolutionary models. 
The BCAH98 models imply a mass of -1.0 Mq (Fig- 
ure rT 2[l. We also compare it to th e evolutionary mod- 
els of lD'Antona k Mazzitelia (119941) whic h yield a mass 
of -0.8 Mq. and iPalla k Stahlerl ([1999D . which yield 
—0.9 Mq. Given these moderate systematic differences, 
we adopt a mass of 0.9±0.1 Mq. Note that this is sig- 
nificantly higher than the value of 0.6±0.1 Mq that has 
been cited in previous work using the later Ml spectral 



type. 

4. DISCUSSION 

4.1. Did GSC 06214-00210 b Experience a Scattering 

Event? 

As described in EJTJ several explanations can account 
for the observed population of planetary-mass compan- 
ions orbiting stars at several hundred AU. Unfortunately, 
these models make few unique or testable predictions, 
and there is considerable debate about formation scenar- 
ios in the literature. Here we focus on the viability of 
one explanation, planet-planet scattering, in the context 
of our results for GSC 06214-00210 b. 

The formation of closely packed planetary systems 
with two or more giant planets naturally leads to dy- 
namical interactions which can significantly alter the 
orbits of one or more components. These close en- 
counters produce a wide range of outcomes including 
(but not limited to) collisions between planets, accre- 
tion onto the star, rearrangement to new stable or quasi- 
stable orbits, scattering to close separations, scattering 
to highly eccentric wide orbits, and complete ejection 
("ionization"). There is growing observational evidence 
that scattering is an important phenomenon in extra- 
solar planetary systems. It has been invoked to ex- 
plain the population of planets wit h high eccentricities 
found in radial velo c ity surveys (e.g..lRasio k FordllT996l : 
iFord k Rasiol 120081 : Uuric k Tremaine!l2008t ). the distri- 
bution of Rossiter-McLau ghlin spin-orbit measurements 
(|Morton k Johnsonll2011l) . and, more recently, the pos- 
sible abundance of unbound gas gian ts inferred from mi- 
crolensing surveys (jSumi 

We examine a hypothetical scattering event for 
GSC 06214-00210 b by assuming it was formed much 
closer than its present location through conventional 
means (core accretion or disk instability) and was ejected 
to a large orbit through the gravitational interaction with 
another massive body. Monte Carlo simulations of scat- 
tering events show that, given unequal planet mass ra- 
tios, the lower-mass planets show a much stronger prefer- 
ence for o utward scattering than t h e more massive com- 
ponents (|Veras k Armitagd [200l IFord k Rasiol 120081) 
with little dependency on which planet initially had the 
wider orbit. If GSC 06214-00210 b was ejected to its 
present location as a result of such an event then it is 
likely that another object at least as massive (—14 Mj up ) 
also formed in that system. 

Ejections of planets onto wide orbits are not un- 
common outcomes in simulations so the discovery of a 
planetary-mass object at —300 AU is perhaps not sur- 
prising. The presence of a disk around GSC 06214- 
00210 b, however, implies that if it underwent scattering 
then its disk was not destroyed during point of closest 
approach with the scatterer. To our knowledge there are 
no studies that specifically investigate the survivability 
of circumplanetary disks during scattering events. Here 
we qualitatively examine the two relevant length scales — 
the separation at closest approach during a scattering 
event and the circumplanetary disk radius — to assess 
the likelihood of a disk surviving such an encounter. 

The criterion for orbital stability betwee n two co- 
planar planets on circular orbits was found by iGladmanl 
(p9l to be 
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(1) 



where A—a2 — a\ is the difference between the initial 
scmimajor axes of the planets and Rh,m is the mutual 
Hill radius (jMarchal fc Bozislll982l ): 



-Rh.m = 



mi + m 2 
3M* 



1/3 



ai + a 2 



(2) 



This radius defines the region in which the gravitational 
force between two bodies is larger than the force on them 
due to the star. The critical range for stability A cr can 
be divided into an outer region encompassing weak in- 
teractions, which can result in simple rearrangements 
of the system architectures with ne w quasi-stable orbits 
(|Fordll2001l : IVeras fc Armita gc 2004), and an inner region 
where strong interactions occur, which is characterized 
by chaotic events made up of collisions and ejections. 
This boundary must be found empiri cally from Monte 
Carlo simulations (see [Gladmanj 119931 for the empirical 
relation for equal mass planets). 

To test the critic al range for stability A cr , 
IChambers et all (1996) ran simulations of two in- 
teracting planets with various separations and confirmed 
that their orbits were stable for at least 10 7 yr when 
A>2v / 3 i?H,M, and when A<2y / 3 i?H,M they eventu- 
ally experienced a close encounter with a distance at 
closest approach of <1 Rh,m- Although the planet 
masses in that study were smaller than mass scales 
for a GSC 06214-00210-lik e system a simi lar study for 
~Jovian mass planets by iMarzaril (|2002f) resulted in 
similar close approach scales of <1 Rh,m and a large 
frequency of ejections. Most scattering studies focus 
on various outcomes as a function of the initial planet 
spacings without discussing the distances of closest 
approach from the simulations, which is the relevant 
length scale we are interested in here. The exception 
is for collisions, defined by various authors to be the 
outcome when A becomes less than about the sum of 
the planetary radii. Given the relatively large frequ ency 
of collisions in some simulations (e.g., iFordl [20011 find 
that collision rates can reach tens of percent depending 
on the planet radii and semimajor axes), the distances 
of closest approach appear to be quite small, which is 
consistent with the se para tion of <1 -Rh .m found by 
IChambers et all (|1996f) and IMarzaril (|2002l )l 

Circumplanetary disks are thought to form around 
young gas giant planets from accretion of circumstellar 
disk material once a planet has opened a gap in the proto- 
plane tary disk (e.g., iLubow et al.l 119991 ; iWard fc Canupl 
2010). In our Solar System the regular satellites of 
the giant planet s are fossil records of these structures 
(jCanup k, Wardl l2002f l and the excess optical emis- 
sion from Fomalhaut b has been attributed to scat- 
tering f rom an eno r mous disk around that extrasolar 
planet (|Kalas et al.l I2008D . The outer radii of these 
disks have been investigated by several authors using 
analytical argume n ts and hydrodynamica l simul ations. 
Guillen fc TrUlm! (fl998lV lAvliffe fc Bate! (pOffflh . and 
Martin fc LubowTd^Oll') obtain similar results of ~0.3- 



0.4 i?H, which roughly corresponds the centrifugal radius. 
(Note that the individual Hill radius Rh is the limiting 
case of the mutual Hill radius Rh,m when m,2 goes to 



zero and a\ =0,2- If to 2 = c m\ and a\ ~ a 2 then i?H,M 
and i?H differ by a factor of (1 + c) 1 / 3 .) 

Since Hill radii depend on the planet masses and semi- 
major axes, we consider the 1:1, 1:2, and 1:5 cases of 
planet mass ratios to better understand the relative mag- 
nitudes of the scales of disk disruption and scattering 
closest approach. For an equal-mass coplanar scattering 
event (to2=toi), the maximum separation at which the 
disks will not be disrupted is (rdi s k i+^disk 2)~0.7 Rh- 
Assuming 01^02 during a close encounter, the conversion 
from mutual Hill radius to the point-particle Hill radius 
is a factor of 1.26, so the approximate scale of closest 
approach is <1.26 Rh- For the case where m2=2mi, the 
disk disruption scale is ~ 0.8 Rh and the distance of clos- 
est approach is <1.4 i?n- Likewise, when m2=5mi, the 
disruption scale becomes ^0.95 i?n and the distance of 
closest approach is <1.8 Rh- 

These results suggest that the disruption of circum- 
planetary disks may be common in scattering events, es- 
pecially if the distance to closest approach is less than 
about one half of a mutual Hill radius. The degree 
to which a disk is affected by these interactions (ie, 
disk truncation vs. complete destruction) likely depends 
on the disk sizes, planet mass ratio, relative disk in- 
clinations, relative planet speeds, and impact param- 
eter. Nevertheless, we interpret the retention of the 
disk around GSC 06214-00210 b as evidence against a 
past scattering event, although detailed simulations are 
needed. 

Further insight into the formation of wide, low-mass 
companions can be gleaned from the architectures of 
other low-mass ratio systems (see fc|4.2[) . In a two-planet 
scattering event involving an ejection, energy conserva- 
tion requires the remaining bound planet to have a final 
semimajor axis greater than one-hal f of its initial semi- 
major axis fe.g.. lFord fc Rasiol[2008l ). Since none of the 
wide planetary-mass companions discovered so far have 
an observed scatterer present in the system, the scat- 
tering hypothesis appears unlikely, although the dearth 
of published detection limits for additional companions 
in these systems hampers a more quantitative assess- 
ment. The notable exception is deep AO imaging o f 
the 1RXS J1609-2105 system bv lLafreniere et al.l pOlOh . 
and somewhat s hallower AO imaging of the GSC 06214- 
00210 system bv llreland et al.l (|2011li . 

Finally, we note that if GSC 06214-00210 b was scat- 
tered to a large orbit then we can constrain the final 
location of the hypothetical scatte rer using the aper - 
ture masking detection limits from iKraus et al.1 (2008) 
combined with the direct imaging detection limits from 
llreland et all (|2011[ ). As shown in Figure [131 the de- 
tection limits exclude objects with the same mass as 
GSC 06214-00210 b at projected separations >170 AU, 
30 Mj up objects (roughly twice the mass of GSC 06214- 
00210 b) >60 AU, and 40 M Jup objects >5 AU. 

4.2. Accretion in Bound and Free-Floating 
Planetary- Mass Objects 

The discovery of strong accretion in GSC 06214- 
00210 b prompts the broader question of how com- 
mon this phenomenon is in other objects. Although 
the census of young planetary-mass companions has 
been steadily growing over the past several years, the 
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moderate-resolution spectra needed to detect infrared 
emission lines have only been obtained for a fraction of 
these objects. 

Six companions with masses 5,20 M,i un and age s 
<10 Myr are known: SR 12 C (iKuzuhara et al.1 120111) . 
CHXR 73 B dLuhman e tal||2006D, 1 RXS J1609f 2105 b 
(lLafreniere et al.l I2008T). ~D H Tau b (|Itoh et all [2o5§h . 
CT ChaB (ISchmidt et al.l200l . andGSC 06214-00210 b 
(note that the uncertainties in masses and ages can 
be rather large). We have ignored companions to 
young brown dwarfs (e.g., 2M1207 Ab) since the for- 
mation of those systems was probably different from 
that of stars and wide companions. Low-resolution 
spectra have been published for most of these compan- 
ions, but moderate-resolution spe ctra have only been ac- 
quired for 1RXS J16 09-2105 b (lLafreniere et al l 12001 
lLafreniere et~aII[20Tf3 ) . CT Cha B (jSchmidt et aljl2008D . 
and GSC 06214-00210 b (this work). CT Cha B and 
CSC 06214-00210 b both show Pa/3 emission, but there 
is no evidence of accretion or thermal disk emission out 
to 4 fan for 1RXS J1609-2105 b. 

Two out of three very low-mass companions there- 
fore show evidence of a disk. The sample size is 
minute, but this could hint that a large fraction 
of these objects are accreting. Moderate resolution 
near-infrared spectroscopy of the remaining compan- 
ions would help address this matter. Intriguingly, 
the low-mass (~25 A/j up ) brow n dwarf companion 
GQ Lup B also show s Pa/3 emission (ISeifahrt et al.ll2007l : 
ISchmidt et all 120081: iLavigne et al.1 12009ft . which may 
support this notion. 

With a large enough sample size, the presence of accre- 
tion or thermal disk emission from young planetary-mass 
companions might be used to learn about their formation 
mechanisms. If this population formed like free-floating 
planetary-mass objects then they should both share sim- 
ilar physical properties such as disk frequencies. Such a 
result would be a sign of a common formation mecha- 
nism. In the future, a similar analysis might be applied 
to extrasolar planets found at close separations (e.g., 
<100 AU) and those at wide separations (>100 AU) to 
test formation mechanisms, perhaps in the context of 
core accretion and disk instability. 

5. CONCLUSION 

Our J- and ii-band spectroscopy of GSC 06214- 
00210 b reveals it is a late-type object (L0±1) with sev- 
eral signatures of youth, including very strong Pa/3 emis- 
sion. The discov ery of accretion in t his object confirms 
the suggestion bv llreland et al.l (|201lD based on its K—L' 
color that it possesses a circumplanetary disk. The accre- 
tion rate of 10~ 10 7 Mq yr" 1 is higher than other objects 
with comparable masses but is consistent with the scatter 
in accretion rates at larger masses. Atmospheric model 
fits to our spectrum yield relatively warm temperatures 
(2700-3000 K), which disagrees with the cooler predic- 
tions from evolutionary models (2200±100 K). With our 
new luminosity measurement of -3.1±0.1 dex, we refine 
the predicted mass of GSC 06214-00210 b to 14±2 M Jup , 
making it the lowest-mass companion to harbor a disk. 
In addition, our optical and near-infrared spectroscopy of 
the primary indicate an earlier spectral type of K7±0.5 
than previously reported. We revise the mass estimate of 
GSC 06214-00210 to 0.9±0.1 M Q based on our updated 



temperature and luminosity. WISE photometry of the 
primary reveals a marginal (2-a) excess at 22 /im. 

GSC 06214-00210 b is one of only a handful of known 
companions orbiting stars at several hundred AU with 
masses straddling the brown dwarf /planetary- mass limit. 
While we cannot unambiguously distinguish the forma- 
tion mechanism of GSC 06214-00210 b, we suggest that 
planet-planet scattering is an unlikely explanation. The 
small distance at closest approach between giant plan- 
ets in a scattering event intimates that circumplanetary 
disks are probably disrupted during such encounters. If 
this conclusion is bolstered by more detailed dynami- 
cal simulations then another explanation — perhaps in 
situ formation — must be relevant for the population of 
planetary-mass companions on wide orbits. This also im- 
plies that the young free-floating planetary-mass objects 
with disks are probably not scattered planets and instead 
represent the low-mass tail of brown dwarf formation. 
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TABLE 1 
Properties of the Primary 
06214-00210 



GSC 



Property 



Value 



Observed 



SpT 

r' (CMC14) 
i (DENIS) 
J (DENIS) 
K s (DENIS) 
J (2MASS) 
H (2MASS) 
K s (2MASS) 
3.4 iim (WISE) 
4.6 /im (WISE) 
12 (WISE) 
22 (WISE) 



K7±0.5 
0.38 ± 0.07 L Q 
11.94 mag 
11.08 ± 0.03 mag 
10.05 ± 0.06 mag 
9.31 ± 0.06 mag 
9.998 ± 0.027 mag 
9.342 ± 0.024 mag 
9.152 ± 0.021 mag 
9.094 ± 0.024 mag 
9.111 ± 0.021 mag 
8.985 ± 0.035 mag 
8.315 ± 0.271 mag 



Estimated 



Distance 
Age 
Mass 
T e ff 



145 ± 14 pc 
5 ±2 Myr 
0.9 ± 0.1 Mq 
4200 ± 150 K 



TABLE 2 

Properties of the Companion GSC 
06214-00210 B 



Property 



Value 



Observed 



SpT 


L0±1 


Proj. Sep. 


2'.' 2 (320 ± 30 AU) 


EW(Pa.0) 


-11.4 ± 0.3 A 


log(£Bol/£©) 


-3.1 ± 0.1 


J (MKO) a 


16.25 ± 0.04 mag 


H (MKO) a 


15.55 ± 0.04 mag 


K (MKO) a ' b 


14.94 ± 0.03 mag 


L' a 


13.75 ± 0.07 mag 



Estimated 



Mass (Evol.) 
T eS (Evol.) 
T eff (BT-Settl) 
log(L(Pa/3)/L ) 

log(Lacc/i©) 

log(M/Af yr" 1 ) 



14 ± 2 M Jup 
2200 ± 100 K 
2700 ± 200 K 
-6.14 ± 0.08 
-4.4 ± 1.3 
-10.7 ± 1.3 



a Photometry from llreland et all J2011I1 . 

b Weighted mean a nd uncertainty of fo ur 
measurements from Ireland ct al. (2011). 
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Fig. 1. — Mcdian-collapscd J-band OSIRIS cube showing GSC 06214-00210 and its companion. The primary is marked with a white 
star. The spaxel scale is / /05 pix -1 and the separation between the primary and the companion is 2'.'2. 
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Fig. 2.— Flux-calibrated J and H band spectra of GSC 06214-00210 b. The strong emission line at 1.282 ^um (inset) is Pa/3 (EW=- 
11.4±0.3 A). The spectrum exhibits absorption features typical of late-M/early-L spectral types including FeH, K I, and H2O. Spectral 
measurement uncertainties are shown at the bottom in gray. 
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Fig. 3. — GSC 06214-00210 b (blue) compared to M8— L2 field objects (left panel) and members of Upper Scorpius (right panels). 
GSC 06214-00210 b resembles the L0.5 field object, although the depth of the K I lines are noticeably weaker in the younger object. The J 
band of GSC 06214-00210 b is well matched to M9 — L2 USco objects an d the H band to M9-L1 members; we adopt a spectral type of L0±1. 
Field objects are from the IRTF Spectral Library HHavner e t al. 2009). From top to bottom they are Gl 752 B (M8), DENIS-P J1048-395 6 
(M9), 2MASS J0746+2000AB (L0.5), 2MASS J0208+2542 (LI), and Kelu-1 AB (L2). The USco spectra are from ILodieu et ail (I2008T) . 
From top to bottom they are USco J155419-213543 (M8), USco J160847-223547 (M9), USco J160737-224247 (L0), USco J160723-221102 
(LI), and USco J160603-221930 (L2). In the left panel the spectrum of GSC 06214-00210 b was smoothed to match the resolving power 
of the field objects (i?~2000). In the right panels the spectrum of GSC 06214-00210 b and the USco spectra were smoothed to i?,~1000. 
The spectra are normalized to 1.29-1.32 fim in the left and middle panels and to 1.60-1.65 fim in the right-most panel. 
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Fig. 4. — BT-Settl-2010 atmospheric model fits (red) to our flux calibrated OSIRIS spectrum (black). J (top), H (middle), and joint J+H 
(bottom) segments are fit separately; the best-fitting models are labeled. Contour plots map the average reduced x 2 values and inferred 
radii for each spectral segment based on Monte Carlo simulations of the data, which account for spectral measurement and photometric 
errors. We exclude the Pa/3 emission line in the fits (gray shaded region). Surface gravity is poorly constrained compared to effective 
temperature. The inferred temperature of 2700-3000 K is significantly higher than evolutionary- model predictions of ~2200 K (Figure |5). 
The reduced x 2 contours represent values of 8, 10, 15, and 20, and the radius contours represent values of 1.1, 1.2, 1.5, and 1.8 i?j up . The 
best-fit model is plotted as a yellow star. Note that the best-fitting model to the joint J+H spectrum is at the edge of the grid. 
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Fig. 5. — Evolutionary m odel-predicted tempe ratur e (top) and ma s s (bot tom) of GSC 06214-00210 b. The l eft panels shows th e Lyon 
Cond and Dusty models of Baraffc ct al. (2003) and Chabricr ct al. (2000), while the right panels shows the Saumon & Marley (2008) 
clear and cloudy (/ so( j=2) variants. In the top panels the models are interpolated onto a grid of constant temperatures. All models yield 
an effective temperature of ~2200±100 K for GSC 06214-00210 b with little variation between the cloudy and clear prescriptions or the 
different models. The Lyon models and the Saumon & Marley models predict masses within 1 Mj up from each other, with the cloudy 
versions yielding slightly higher masses than the clear variants. We adopt a mass of 14±2 Afj up for GSC 06214-00210 b. 
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Fig. 6. — Mass accretion rates (M) versus stellar mass. The data are fr om Natta et al. (2006, accret i on ra tes derived using Pa/3 lines are 
shown in blue while those using Br7 are in orange). IHerczee et al, (2009, green), and Muzerolle et al. ( 2005I. red). Open inverted triangles 
represent upper limits. Although the various samples represent different ages and methodologies used to derive mass accretion rates, the 
data clearly show the strong dependency of accretion rate on mass. The accretion rate of GSC 06214-00210 b appears to be somewhat 
higher than for low-mass brown d warfs. The uncerta inties shown here in the accretion rate for GSC 06214-00210 b ignore errors in the 
conversion from Lp a ^ to L a cc from Natta et al. ( 2004) since we are making a direct comparison to their sample. The uncertainty increases 
to ±1.3 dex when the errors are included. 
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Fig. 7.— Comparison of the Pa/3 B1V (top) and line luminosity (bottom) for GSC 02614-00210 b to the INatta et all pOOBT ) sample. 
Filled circles represent detections and open inverted triangles represent upper limits. The position of GSC 02614-00210 b is marked with 
a yellow star. 



20 




Fig. 8. — SXD spectrum of the primary GSC 06214-00210. No Pa/3 or Br7 emission lines are evident (upper right insets), which is 
consistent with mid-infrared photometry showing no indication of a disk. The bottom left inset shows that the SpeX/prism low resolution 
spectrum (red) agrees well with the SpeX/SXD moderate resolution data. 
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Fig. 9. — Left: UH 2.2m/SNIFS spectrum of the primary (red) compared to dwarf stars from the Pickles (1998) spectral library. The 
SNIFS spectrum is a composite of blue (3300-4900 A) and red (5170-8700 A) channels. GSC 06214-00210 is we ll matched by the K7 
template. Right: SXD spectrum (red) compared to dwarf stars from the IRTF Spectral Library (Rayncr ct al. 2009). From top to bottom 
the objects are HD 45977 (K4), HD 36003 (K5), HD 201092 (K7), HD 19305 (M0), and HD 42581 (Ml). K4-K7 objects are good matches 
to our near-infrared spectrum of the primary. We adopt a spectral type of K7±0.5 for GSC 06214-00210. 
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Fig. 10. — Spectral energy distributions of GSC 06214-00210 and its companion. The T e fr=4200 K/log g=i.O Phoenix-Gaia model match 
the 0.6-12 /im photometry of the primary. For the companion the T of f=2700 K/log 3=4.0 and T eff =2200 K/log g=4.0 BT-Settl-2010 
models are plotted. The warmer temperature is from fitting atmospheric models and the cooler temperature is the evolutionary model 
prediction. Both models are flux calibrated to the J band photometry of the companion. The excess flux in L' , and perhaps also in K are 
likely caused by thermal emission from a circumplanetary disk. Uncertainties in the photometry are smaller than the symbol sizes except 
for the WISE 22 /an point. The slight excess at 22 fim seen in the primary may be caused by a circumstellar disk. 
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Fig. 11. — W1-W4 vs J-Ks diagram for field K7 stars (black filled circles) and Upper Scorpius members with spectral types of G, K, 
and M (filled red circles). Field stars have WI—W4 colors near 0.0 mag, while many of the redder USco members show an excess at these 
wavelengths. WISE photometry of GSC 06214-00210 (blue star) suggests an excess at the 2-<r significance level. 
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Fig. 12. — HR diagram for GSC 06214-00210 and its companion. The evolutionary models of lBaraffe et alj l|1998l ) with L mix =H p (black) 
are overplotted with 1, 5, 10, and 100 Myr isochrones (dashed lines) and iso-mass tracks (solid lines) from 1.4-0.02 Mq. The 1 Gyr main 
sequence isochrone is shown as a thick line and the hydrogen burning minimum mass (HBMM) of ~0.072 Mq is shown as a dotted line. 
The L m ; x =1.9iJp isochrones for 1 Myr, 5 Myr, 10 Myr, and 1 Gyr are shown in green. The position of the primary indicates a mass of 
~1 Mq; our adopted mass of 0.9±0.1 Mq is based on an average from several evolutionary models (see text for details). The position of 
GSC 06214-00210 b (red) is shown for the warmer atmospheric model-inferred temperature (2700 K) and the cooler evolutionary model- 
derived temperature (2200 K). The warmer temp erature is inconsistent with the young age of the system. Gray circles show members of 
Upper Scorpius from Preibisch & Mamaick (2008). 



25 



60 
50 



Q- 



to 
to 
03 











Aperture Masking 
Direct Imaging 




GSC 0621 4-0021 Ob 



_i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i_ 



100 200 300 

Projected Separation (AU) 



Fig. 13. — Constraint s on additiona l com panions around GSC 06214-00210 based on aperture masking results from lKraus et aLl (2008) 
and direct imaging from Ireland ct al. (2011). Objects with masses twice that of GSC 06214-00210 b are excluded at projected separations 
>60 AU. The aperture masking limits exclude brown dwarfs with masses >35 Mj up at projected separations >5 AU. The uncertainty in 
the projected separation of GSC 06214-00210 b (30 AU) is dominated by the distance error. 



